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Theoretical study of 7 Doradus pulsations in pre-main sequence stars 
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The question of the existence of pre-main sequence (PMS) 7 Doradus (7 Dor) has been raised by the observations of 
young clusters such as NGC 884 hosting 7 Dor members. We have explored the properties of 7 Dor type pulsations in a 
grid of PMS models covering the mass range I.2M0 < A/, < 2.5Mq and we derive the theoretical instability strip (IS) 
for the PMS 7 Dor pulsators. We explore the possibility of distinguishing between PMS and MS 7 Dor by the behaviour 
of the period spacing of their high order gravity-modes (g-modes). 
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1 Introduction 

The variability of 7 Dor stars was identified as due to pul- 
sations by Balona, Krisciunas & Cousins (1994), and the 
features of the stars belonging to this new class of pulsators 
were defined by Kaye et al. (1999). They are late A and 
F-type stars covering a part of the Hertzsprung-Russel dia- 
gram (HRD) between 7200 — 7700 K on the zero-age main 
sequence (ZAMS) and 6900 - 7500 K above it (Handler 
1999), between the solar-like stars and the S Scuti (S Set) IS. 
They are located between stars with a deep convective enve- 
lope (CE) and stars with a radiative envelope, in the region 
of the HRD where the depth of the CE changes rapidly with 
the effective temperature of the star. 7 Dor pulse with high 
order 5-modes in a range of periods between 0.3-3 days. 
The excitation mechanism proposed by Guzik et al. (2000) 
using the frozen-convection assumption is a modulation of 
the radiative flux at the base of the CE. This mechanism 
was revisited by Dupret et al. (2005) using a time-dependant 
convection (TDC) treatment. 

From an observational point of view, the limits of the 
7 Doradus IS have been lastly established by Handler & 
Shobbrook (2002) (HS02 hereafter) and in the rest of the 
paper these Umits will be adopted to define the 7 Dor IS. 
Since the depth of the CE plays a major role in the driv- 
ing mechanism of 7 Dor pulsations, the theoretical predic- 
tions of stability are very sensitive to the parameter ajv/LT 
defining the travel length of convective elements in the clas- 
sical mixing-lenght treatment (MLT) of convection (Bohm- 
Vitense 1958). 

Number of observational efforts have been devoted to 
the search of PMS 7 Dor pulsators. Saesen et al. (2010) 
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Fig. 1 PMS and MS (resp. full & dashed thin lines) evolu- 
tionary tracks. Thick straight lines correspond to the limits 
of the observational 7 Dor IS (HS02). Grey points repre- 
sent the bona fide 7 Dor stars from Henry, Fekel & Henry 
(2005). 



found 6 multiperiodic A and F-type stars with a mean fre- 
quency between 0.2 d~^ and 3 d~^ during their multisite 
observation campain on the young open cluster NGC 884 
(age - 12.8 Myr - Slesnick et al. 2002). Zwintz et al. (2009) 
searched for PMS pulsators in a young open cluster (NGC 2264,| 
age - 3-10 Myr - Sung et al. 2004; Sagar et al. 1986) but did 
not find any 7 Dor which were confirmed cluster members. 

The above mentioned theoretical works on 7 Dor stars 
systematically studied MS models stability. However, as shown| 
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in Fig.[T] MS and PMSU evolutionary tracks cross the obser- 
vational 7 Dor IS. The presence of stars at different evolu- 
tionary phases in this region of the HRD raises some ques- 
tions: since the internal structure of PMS stars is differ- 
ent, can we expect 7 Dor pulsations in such stars? Could 
these pulsations be used to distinguish between PMS and 
MS 7 Dor? 

To answer these questions, we performed an adiabatic 
and a non-adiabatic asteroseismic analysis on a grid of PMS 
and MS models between 1.2 and 2.5Mq computed with 
the stellar evolution code CLES (Scuflaire et al. 2008a). 
The adiabatic and non-adiabatic computations have been 
done respectively with the LOSC (Scuflaire et al. 2008b) 
and MAD (Dupret 2001) codes. 



2 Internal structure and adiabatic study 
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Fig. 2 Top panel: Evolution of a 1.8Af0 star in the HRD 
(PMS phase in full black and MS phase in dashed grey). 
Middle panel: Variation of its CC mass from the PMS phase 
to the early MS phase (Xc — 0.68). Bottom panel: Evo- 
lution of the £ — 1 modes period spacing as a function of 
the effective temperature of the star from the PMS to the 
Terminal Age Main Sequence (TAMS). 
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Fig. 3 Evolutionary tracks of different models crossing 
the observational 7 Dor IS (thick grey lines - HS02). Left 
panel: I.SMq evolutionary track whose PMS (full line) and 
MS (dashed line) phases intersect inside the IS (circle: PMS 
model - triangle: MS model). Right panel: Evolutionary 
tracks for 1.9 and 2.1A^q showing the same HRD location 
for a 2.IM0 PMS model (circle) and a model of 1.9Mq at 
the end of its MS phase (triangle). 



A star approaching the MS from the Hayashi track has 
already a radiative core that continues to contract. The in- 
crease of density due to this contraction leads to the increase 
of the central temperature (Tc). This phenomenon continues 
until Tc is high enough (^ 1.7 10^ K) to start the nuclear re- 
actions of the CN subcycle. Because of the high dependence 
on temperature of the ^^C(p, 7)"N(;3+iy)i^C(p, 7)^'*N nu- 
clear reaction rate (oa T^^) a convective core (CC) appears. 
The fraction of mass of this CC changes as the star evolves 
toward the ZAMS (Fig. |2] - middle panel) and for typical 
7 Dor stellar masses, the CC remains during the MS. The 
onset of the CN subcycle appears in the PMS evolutionary 
track as a kind of loop with a minimum of luminosity after 
a first maximum (Fig.|2]- top panel). 

As shown in Fig. [T] the phase at which the star crosses 
the IS during the PMS changes with the stellar mass. While 
low mass models have already developed a CC, more mas- 
sive models are still contracting with a fully radiative core. 

The properties of the g-modes spectrum is determined 
by the matter stratification in the star, which is described by 
the Brunt-Vaisala frequency TV: 



N'^g 



1 dlnP dlnp 



Fi dr 



dr 



(1) 



' we consider as PMS models those before the onset of the stationary 
central H-burning 



with g the local gravity, p the local density, P the local pres- 
sure, r the local radius and Fi the first adiabatic exponent. 
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From the first order asymptotic theory (Tassoul 1980), 
the period of a g-mode with a radial order k and a degree £ 
in a star with a CC and a CE is given by: 



Pk = 



vw^)i:: ^ 



dr 



{2k + 1) 



(2) 



with ri and r2 the limits of the g-mode cavity defined by 



cri < N' 



q2 



where Ug is the g-mode frequency and Se 



the Lamb frequency for modes of degree £ (Fig. |4|. The de- 
pendence of g-mode periods on the behaviour of N in the 
central region of the star allows us to expect a clear differ- 
ence between the seismic properties of PMS and MS 7 Dor. 

2.1 Comparison of the models internal structures 

The PMS evolutionary track of a I.SMq star crosses at dif- 
ferent points the corresponding MS track. In particular, both 
tracks cross in the middle of the 7 Dor IS (Fig. [3] - left 
panel). By comparing these two models at the same loca- 
tion in the HRD we eliminate undesirable effects on the stel- 
lar structure coming from different effective temperatures or 
luminosities. The internal structures of these PMS and MS 
models are shown in Fig. |4] where we plotted their N and 
Se=i frequencies as a function of the relative radius. 

Because of the same radius and mass, PMS and MS 
models show similar behaviour of N in the outer layers and 
the bases of the CE are located at the same depth (Fig. |4]- 
left panel). Both N profiles present also a bump in the inner 
layers, due to the density distribution. However, the central 
layers of the two models are very different. The PMS model 
has only a small CC {N^ < 0) due to the onset of the CN 
subcycle while the CC of the MS model is larger The main 
difference between the two N profiles lies in the sharp fea- 
ture located at the limit of the MS convective core. This 
peak is due to the presence of an important mean molecular 
weight gradient (V^) at the limit of a receding CC. 

As mentioned above, more massive PMS models cross 
the 7 Dor IS during an earlier phase than lower mass ones. 
We considered models with CC overshooting in order to 
have MS models with high enough luminosity and low enoug 
effective temperature crossing a more massive PMS track 
(Fig. |3]- right panel). A 2.1^/© quasi chemically homoge- 
neous PMS model that still has a radiative core has the same 
HRD location as a 1.9Mq evolved MS model. The differ- 
ent masses but same radius of the two models lead to differ- 
ences in their density profiles that determines the behaviour 
of N in the envelope, and therefore the different depth of 
their CE (Fig.|4]- right panel). 

2.2 Adiabatic study - Period spacing 

From Eq. 2, the period spacing between two g-modes with 
consecutive radial orders and same degree can be written as: 

27r2 
AP = P,.+i - P, = ————^^^ (3) 
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Fig. 4 Propagation diagram. PMS (black) and MS (grey) 
N and "£ = 1"-Lamb (dotted lines) frequencies for mod- 
els with the same mass (left) and for models with different 
masses (right). The straight thick lines represent the propa- 
gation of g-modes having the same radial order in the PMS 
(black) and MS (grey) models. 

Fig. I2] (bottom panel) represents the evolution of the 
period spacing from the PMS to the TAMS for a l.SM© 
star The period spacing presents a clear variation which is 
strongly dependant on the evolution of the CC. 

For models with the same effective temperature, the max-| 
imum difference between PMS and MS period spacing is 
around 2000 seconds and corresponds to the maximum mass| 
fraction of the CC (Fig. |2] - point E). The period spacing 
value could be, in principle, used as a first discriminant be- 
tween PMS and MS 7 Dor. However the period spacing of 
PMS stars can be of the same order as the MS 7 Dor one 
(Fig. ID - bottom panel - points A & C) and we cannot al- 
ways use this value to determine if a 7 Dor star is in its 
PMS or its MS evolutionary state. 

Nevertheless, a sharp variation of N such as the one due 
to the Vp at the border of the CC let a clear asteroseis- 
mic signature: the oscillation of the period spacing around 
its mean value (Miglio et al. 2008 and references therein), 
'while MS models can present that V^^, PMS ones are al- 
most chemically homogeneous, i.e. their N profile is quite 
smooth. Therefore the PMS period spacing does not clearly 
change with the radial order (Fig.|5]l. 

3 Non-adiabatic study - Stability analysis 

Because the adiabatic study does not take into account the 
excitation and damping of modes, we performed a non-adiabatic| 
analysis. The theoretical IS were calculated for £ = \ and 
£ ^ 2 modes with the non-adiabatic code MAD includ- 
ing TDC treatment (Grigahcene et al. 2005) on a grid of 
1.2 — 2.5M0 stellar models computed with the following 
phyical inputs: 

- OPAL2001 equation of state (Rogers & Nayfonov 2002) 
and OP opacity tables (Badnell et al. 2005) completed 
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Fig. 5 Period spacing structure for £ ^ 1 modes as a func- 
tion of the radial order of the modes for PMS (black trian- 
gles) and MS (grey squares) models with the same mass 
(M = 1.8Af0) presented in subsection 2.1 and Fig.|3] 




Fig. 6 7 Dor theoretical IS for MS models with over- 
shooting (aov — 0.20, middle grey) and without (dark grey) 
and for PMS models (light grey). Straight lines correspond 
to the observational 7 Dor IS (HS02) and the thin dotted 
lines are the evolutionary tracks for models between 1.4 and 
2.5Mq. 



at low temperature (T < lO^K) by the tables provided 
by Ferguson et al. (2005). The atmosphere tables were 
derived from Kurucz atmosphere models starting at the 
photosphere. We used the abundances from Asplund, 
Grevesse & Sauval (2005) for an initial metal mass frac- 
tion Z = 0.02 and an initial central hydrogen fraction 
X = 0.70. As we used the same non-adiabatic tools as 
Dupret et al. (2004), we chose the MLT parameter value 
they needed to tune their theoretical MS 7 Dor IS to the 
observational one, i.e. umlt — 2.00. They showed that 
a smaller value of the MLT parameter would shift the 
whole theoretical IS to lower effective temperatures. Fi- 
nally, we computed models with and without convective 
core overshooting {aov — 0.20; 0.00). 



Our first results are presented in Fig. |6] showing the the- 
oretical IS of high order g-modes for MS and PMS models. 
The location of PMS 7 Dor theoretical IS matches up with 
the MS one. 

4 Conclusion 

We carried out an adiabatic and a non-adiabatic studies on 
a grid of PMS and MS models in the mass range 1.2Af0 < 
Af* < 2.5A^0 to differenciate the asteroseismic behaviour 
of PMS 7 Dor from that of MS ones. 

We pointed out the theoretical existence of PMS high 
order g-mode pulsators in the region of the observational 
7 Dor IS. The theoretical PMS IS has the same edges than 
the MS one and presents a good agreement with the obser- 
vational 7 Dor IS for umlt = 2.00. The only difference 
between MS and PMS IS lies in the fact that even with a CC 
overshooting no MS evolutionary track can reach the upper 
region of the HRD, where unevolved massive PMS 7 Dor 
exist. 

The measurement of the period spacing allows us to 
make the distinction between PMS and MS models by two 
different ways: 

- At fixed stellar parameters, the difference between cen- 
tral internal structures may lead to a significant differ- 
ence between the values of the period spacing. We are 
aware that we should investigate if we can still distin- 
guish PMS stars once we consider stellar parameters un- 
certainties but it is not the purpose of the present paper. 

- the behaviour of the period spacing is also different be- 
tween MS and PMS models. While during the MS the 
important V^ and the evolution of the CC leads to an 
oscillation of the period spacing, the lack of such a V^ 
during the PMS phase leads to a period spacing inde- 
pendant of the mode radial order. 

In a forthcoming paper the difference between PMS and 
MS non-adiabatic g-modes frequency spectra will be pre- 
sented. 
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